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Abstract 

We present a first-principles study of the electronic structures and properties of ideal (atomically 
sharp) LaAlOa/ SrTiOs (001) heterointerfaces and their variants such as a new class of quantum 
well systems. We demonstrate the insulating-to-metallic transition as a function of the LaA103 
film thickness in these systems. After the phase transition, we find that conduction electrons are 
bound to the n-type interface while holes diffuse away from the p-type interface, and we explain 
this asymmetry in terms of a large hopping matrix element that is unique to the n-type interface. 
We build a tight-binding model based on these hopping matrix elements to illustrate how the con- 
duction electron gas is bound to the n-type interface. Based on the 'polar catastrophe' mechanism, 
we propose a new class of quantum wells at which we can manually control the spatial extent of the 
conduction electron gas. In addition, we develop a continuous model to unify the LaAlOs/SrTiOa 
interfaces and quantum wells and predict the thickness dependence of sheet carrier densities of 
these systems. Finally, we study the external field effect on both LaAlOa/SrTiOs interfaces and 
quantum well systems. Our systematic study of the electronic reconstruction of LaAlOa/SrTiOs 
interfaces may serve as a guide to engineering transition metal oxide heterointerfaces. 
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I. INTRODUCTION 



With the advance of techniques to control thin film growth on the atomic scale, the study 
of epitaxial oxide heterostructure, is a rapidly developing area of material .science Q. Due 
to the ability to produce a well-defined single terminated substrate surface [2J, interfaces that 
are nearly atomically sharp can be formed. In many cases, the properties of these interfaces 
turn out to be much richer than those of the bulk constituents. One of the most interesting 
examples is the (001) interface between LaA103 (LAO) and SrTiOs (STO), both constituents 
of which are conventional band insulators in the bulk. Among the intriguing phenomena 
observed at this interface is the presence of a high-mobility quasi two-dimentional electron 



gas [3|, which is sensitive to the LaAlOs film thickness on the nanoscale 



4-61, and can be 
9] 



tuned via an external field . In addition, both magnetism [8[ and superconductivity 
have been observed at LaAlOs/SrTiOs interfaces. 

Due to these interface properties, the LaAlOs/SrTiOs system shows great promise for 



the development of novel applications in nano-scale oxide electronics lOj, 111). However 



though significant efforts have been made in both theory [12Nl8| and experiment Ji3-25|, the 
origin of the new phases emerging at the LaAlOs/SrTiOs interface has not been completely 
resolved yet. In part, this is due to the fact that more than one mechanism may play a 
role in determining the behavior, and furthermore, the dominant mechanism can depend on 
external conditions 

One mechanism that is generally believed to play a role in the formation of the two- 
dimensional electron gas is known as the polar catastophe. This mechanism is a direct 
result of the charge discontinuity that occurs at an abrupt interface between a nonpolar 
(SrTiOs) and a polar (LaAlOs) material. Both materials are perovskite oxides with an 
ABO3 structure that forms alternating planes of AO and BO2 stacked along the (001) di- 
rection. Consequently, there are two types of LaAlOs/SrTiOs interfaces along this direction: 
Ti02/LaO (known as the n-type interface) and SrO/A102 (the p-type interface) [jj, 0]. In 
the ionic limit, SrTiOs is composed of charge neutral atomic layers (SrO)° and (Ti02)°, 
while LaAlOs consists of positively charged (LaO) + and negatively charged (A10)~ atomic 
layers. When LaAlOs is grown epitaxially on the SrTiOs substrate along the (001) direc- 
tion, an internal electric field directed from (LaO) + to (A102)~ is formed through LaAlOs, 
resulting in a potential difference that, mathematically, diverges with increasing LaAlOs 



2 



thickness. In reality, to offset the diverging potential, an electronic reconstruction is ex- 
pected to occur, with electrons transferring from the film surface (or a p-type interface) to 
the n-type interface. The polar catastrophe mechanism (3), [^j has been confirmed in density 
functional theory (DFT) calculations fla . 28-3o|, and could be responsible for the presence 



of the two-dimensional electron gas at the atomically sharp n-type interface (i.e., with no 
defects) . 

However, the polar catastrophe mechanism alone can not explain all of the electronic 
properties of the LaAlOa/SrTiOs interface, such as the observed confinement of the trans- 



3l| . In previous work 28] , 



ferred electrons within several nanometers of the n-type interface 
we have shown that the electrons are bound to the n-type interface as a result of the chem- 
ical environment of the interface, which produces a large hopping matrix element between 
La and Ti. 

In this paper, we use both DFT and phenomenological modeling to provide a more 
detailed and complete picture of the electronic reconstruction at ideal LaAlOs/SrTiOs in- 
terfaces. We predict the behavior of the sheet carrier density as a function of the LaA103 
thickness by constructing a continuous model which approximates the LaA103 as a homoge- 
neous media. We also use both DFT calculations and model calculations to give a simple yet 
quantitative picture of the external field effect that has been experimentally realized at the 
n-type interface Q,0|. In addition, we propose a new class of quantum well systems, based 
on the polar catastrophe mechanism, at which the spatial extent of the two-dimensional 
electron gas can be manually controlled. 

The remainder of the paper is organized as follows. In Section [Til we discuss the technical 
details of our DFT calculations. We present our first-principles results in Section [TTTl Section 
IIII Al gives a brief discussion on the symmetric superlattices. In IIIIB| we briefly describe 
the polar catastrophe at various types of interfaces, and we introduce the new quantum 
well systems. Section IIII CI discusses the computation of the on-site and hopping matrix 
elements and a tight-binding model. An exceptionally large Ti-La hopping is found at the 
n-type interface and its origin is explained. Section IIII Dl studies the behavior of oxygen 
vacancies at the interfaces. Section IIII El examines the thickness dependence of the sheet 
carrier density. Simulations of the external field effect are presented in IIII Fl In Section IIV| 
we present a continuous model which gives a simple yet quantitative description of both the 
sheet carrier density and the external field effect. We conclude in Section IVl 
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II. METHOD DETAILS 



Our calculations are performed using density functional theory within the ab initio su- 



percell plane- wave approach |32|, with the code PWscf in the Quantum-ESPRESSO pack- 
age |52|. We employ the local density approximation (LDA) [33J and ultrasoft pseudopo- 



tentials 



341 ]. which include semi-core states in Sr, Ti and La. The reference configura- 



tions of the pseudopotentials are: Sr 4s 2 4p 6 5s 2 



r v = r * = i.8a ), La 5s 2 5p 8 5d 1 6s 1 - 5 6p - 5 
1.8a ,rP = 1.82a ) and O 2s 2 2p 4 (r s c 



r s c = 2.0a ,rP = 1.8a ), Ti 3s 2 3p 6 3dHs 2 
r s c = r d c = 2.2a ,r p c = 2.0a ), Al Ss^p 1 
1.3ao), where a is the Bohr radius. The 



plane wave basis energy cutoff and charge cutoff are 30 Ry and 180 Ry, respectively. We 
use a Gaussian smearing width of 5 mRy when sampling the Brillouin zone. The fc-grid 
sampling of the Brillouin zone is 10 x 10 x 1 where the z-axis is orthogonal to the interface. 
The convergence of the total energy and total charge density has been checked with fc-grids 
of up to 20 x 20 x 1. Periodic copies are separated by ~15 A of vacuum. We also extend 
the vacuum to 30 A to check the convergence. The effect of the artifical electric fields in 
the vacuum due to the periodic boundary condition, which turns out to be quite small, 
is discussed in Appendix IB II The force convergence threshold is 26 meV/A. In some key 
results we reduce the threshold to 13 meV/A to check the convergence. 

In all our calculations, the interfaces are along the (001) direction so that the z axis is 
perpendicular to the interface. The x and y directions of the simulation cell are subject to 
periodic boundary conditions and their lengths are fixed to the theoretical lattice constant 
of SrTiO"3 a = 3.85 A (1.5% smaller than the experimental value). The atomic coordinates 



are relaxed 



53] until every force component is smaller than the convergence threshold. The 



detailed configurations for different calculations are specified in each section below. 



III. DFT RESULTS 



A. Symmetric double n-type and p-type superlattices 



The simplest approach to studying 
use a symmetric superlattice approach 



;he n-type or p-type SrTiOs/LaAlOa interface is to 



12 



17 



35 



371 ] , as exemplified by Fig. (U which shows 



a simulation cell without vacuum in which a number of SrTiC>3 unit cells are adjacent to 
a number of LaA10"3 unit cells and periodic boundary conditions are imposed. Therefore, 
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FIG. 1: Schematics of symmetric superlattices. a) The double n-type superlattice. b) The double 
p-type superlattice. The interface is highlighted by the dashed line. 

a no n-stoichio metric supercell contains two identical SrTiOs/LaAlOs interfaces: e.g. in 
Fig. [T^l, we have an additional TiC>2 layer in SrTiC>3 and an additional LaO layer in LaA103, 
and in Fig.[T]D, there is an additional SrO layer in SrTiC>3 and an additional AIO2 in LaA103. 
In this way, the superlattice contains two n-type interfaces (TiC^/LaO) or two p-type inter- 
faces (SrO/AlOa). 

The advantage of the superlattice approach is that no vacuum is needed in the simu- 
lation cell, making computation easier. However, due to the imposed symmetry and non- 
stoichiometry of the LaA103 film, this geometry does not result in a polar field, so the 
evolution of the polar catastrophe can not be modeled. In addition, the non-stoichiometry 
of the LaAlOs also imposes a fixed carrier doping: in the ionic limit, an LaO (AIO2) layer 
has a charge of +1 (-1), and there is one extra electron (hole) present in the conduction 
(valence) band of LaA103, which is shared evenly by the two interfaces. Therefore each 
interface is doped by 0.5 electrons (n-type) or 0.5 holes (p-type) per two-dimensional unit 
cell. These values are precisely those needed to fully compensate the polar field of LaA10 3 
[27]. so the symmetric supercell approach is equivalent to studying the properties of the 
interfaces when the LaA103 film is very thick (infinite thickness limit). 

We briefly present results on the symmetric interfaces and highlight some key observations 
and questions that we will answer in later sections. We show in Fig. [2|i the xy-plane 
integrated conduction electron density for the double n-type superlattice, and in Fig [2b 
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FIG. 2: Transferred electron and hole densities integrated over the xy plane. In each panel, the 
integral of the conduction electron/hole density is normalized to unity. The layers are measured 
in units of the SrTiOs lattice constant a. a) The symmetric n-type superlattice. The electrons 
are bound to the n-type interface, b) The symmetric p-type superlattice. The holes reside in both 
the SrTiO"3 and LaAlC>3 layers, c) The n-type interface with 4 u.c. of LaA10"3. The electrons 
decay away from the n-type interface, d) The p-type interface with 5 u.c. of LaA10"3. The holes, 
driven by the polar field through LaA103, diffuse into the SrTiO"3 substrate. Appendix lAl explains 
in detail how to calculate these transferred charge densities. 

the hole density for the double p-type superlattice. The method of how to calculate the 
transferred charge density is detailed in Appendix [S] The integrated densities are 0.49 
electrons per n-type interface and 0.49 holes per p-type interface, showing that the ionic 
picture for LaA103 is highly accurate, and that one should expect 0.5 electrons (holes) at 
the fully compensated n-type (p-type) interface for very thick LaA103 films. 

The density profiles in Fig. [2] display an interesting asymmetry. One can characterize 
the electrons in the double n-type system as being bound to the interface. As the figure 
shows, they are confined in the SrTiC>3 , with the density decreasing away from the interface. 
On the other hand, in the double p-type system, the holes reside in both the LaA103 and 
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SrTiOa, do not show a strong preference for the interface itself, and have a slightly higher 
amplitude in the LaAlOs. 

Our calculated band offsets place the LaA103 valence band maximum ~ 0.1 eV above 
the SrTiOs valence band maximum (p-type) and the LaAlOs conduction band minimum 

n n 

~ 2 eV above the SrTi03 conduction band minimum (n-type) [13|, [15(. A priori, one would 
expect that any introduced carriers would occupy the most energetically favorable band 
edge available, i.e. holes would migrate to the LaAlOs valence band edge and electrons to 
the SrTiOs conduction band edge. Furthermore, the minimization of kinetic energy would 
lead to a relatively uniform distribution of carriers. In Fig. [2J we see that, overall, the holes 
prefer to be in the LaA103 at the p-type interface, while the electrons reside in SrTiOs at 
the n-type interface. For the p-type interface, the small value of the band offset and the 
large hole density result in some "leakage" into the SrTiOs valence band states. 

However, the binding of the electrons at the n-type interface cannot be explained by the 
above arguments. As no polar fields are present in this system to create a potential that 
will bind the electrons, what mechanism overcomes the kinetic cost inherent in localizing 
them to the interface? Whatever mechanism is present is absent for the holes, whose spatial 
profile is relatively uniform in both materials. We provide the answer to this fundamental 
question in Section IIII C\ when we discuss the tight-binding hopping matrix elements across 
both interfaces. 

Before we move onto the next section, we would like to give a short discussion on the 



valence band offset (VBO) of LaAlOs /SrTiOs inter 



aces. The VBO of the n-type interface 



has been calculated by several groups [lOj, |13|, [15|, |37|, |38|. Although all of the computed 
values place the conduction band edge of LaA10 3 significantly higher than that of SrTi0 3 
at the n-type interface, the large range of computed VBO values (shown in Table I), is 
puzzling. To investigate this issue, we use two different approaches to calculate the VBO of 
the n-type interface and the p-type interface. The valence band offset is defined by: 

-&VBO — — -^V K L ) 

where E^ AO and E§ TO are the valence band edges of LaAlOs and SrTiOs, respectively. 
The simulation cell is a symmetric nonstoichiometric superlattice with 12.5 layers of SrTiOs 
and 4.5 layers of LaAlOs. The first method to determine the VBO is to analyze the local 



density of states, referred to as the LDOS approach 
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391 ] . The other method is to use 



the macroscopic average potential, denoted further as the MA approach [15|, |4Cj, |4l| . We 
compare our results with other theoretical calculations and the recent experiment in Table 

m 



TABLE I: The valence band offset (VBO) of LaAlOa/SrTiOa n-type and p-type interfaces in eV. 
Each value is followed by the method used to determine it in parentheses. 





n-type (Ti0 2 /LaO) 


p-type (SrO/A10 2 ) 


present study 


-0.47 (MA) -0.39 (LDOS) 


-0.02 (MA) 0.11 (LDOS) 


Ref. [15] 


0.51 (MA) 0.39 (LDOS) 


0.19 (MA) 0.10 (LDOS) 


Ref. [10] 


1.1 


NA 


Ref. [13] 


-0.15 (LDOS) 


NA 


Ref. [37] 


-0.9 (LDOS) 


NA 


Ref. [38] 


-0.1 (MA) 


NA 


Expt. [42] 


-0.35±0.18 


NA 



We can see from Table [I] that for the n-type interface not only the magnitude of VBO 
differs, but the sign is not unanimous. One possible origin for the range of theoretical values 
is likely that the lattice constants have minute differences (due to different pseudopotentials), 
which causes a variation in strain and probably affects the alignment of valence band edges. 
Another possibility is that we, as well as a few other groups 13|, |l5|, |37|, [38J], use symmetric 
nonstoichiometric superlattices, in which the macroscopic average potentials are flat in both 
LaAlOs and SrTiOs, to perform the calculations. However, if a stoichiometric LaAlOs slab 
is employed [lo| so that a net internal electric field results, the determined VBO will likely 
be different from the symmetric nonstoichiometric case. Although our results for the n-type 
interface agree well with the available experiments, the value as well as variations of the 
VBO for the n-type interface is not a closed subject and needs further work. In contrast, 
the computed values of the VBO for the p-type interface are in general agreement and ~ 
0.1 eV. 



S 



B. The polar catastrophe 



1. n-type, p-type and np-type interfaces 

As discussed in the Introduction, there are two types of LaAlOa/SrTiOa interfaces: the n- 
type (Ti02/LaO) interface and the p-type (SrO/AlC^) interface. If a stoichiometric LaA103 
layer on an SrTiOa substrate also has a SrTiOs capping layer, then we have both n-type 
and p-type interfaces in the same system. Below, we refer to this type of geometry as the 
np-type interface. Experimentally, the n-type, p-type and np-type interfaces have all been 
fabricated. In this study, we simulate all three types of interfaces, the configurations of which 
are illustrated in Fig. [3J For the np-type interface, we use 5 unit cells of SrTi0 3 to simulate 
the substrate and another 5 unit cells of SrTi03 as a capping layer. For both the n-type 
and p-type interfaces, we use 11 unit cells of SrTi0 3 to simulate the substrate. In all cases, 
the thickness of LaA103 is varied from 1 to 7 unit cells. In addition, the first two SrTiOa 
unit cells in the substrate are fixed at the ideal perovskite positions to simulate the bulk-like 
interior of the substrate 28j. The termination of the SrTiOa substrate and capping layers 
is always SrO. As illustrated in our previous work 54j, the TiC>2 termination has a surface 
O-p state which is ~ 0.5 eV higher than the valence band edge of bulk SrTiC>3, while the 
SrO termination does not have this surface state. Since we are interested in the evolution 
of polar fields of LaA103 on very thick SrTi03 substrates, this pure surface state should be 
avoided when simulating the SrTi03 substrate and thick SrTi03 capping layers. However, 
for thin SrTi03 capping layers, the presence of this O-p state on the TiOo-terminated surface 



will lower the critical separation, which has recently been discussed [43 1. 

We observe that an insulating-to-metallic transition occurs in our simulations when the 
LaA103 film reaches a critical thickness, the value of which depends on the system geometry. 
The phase transition can be explained in terms of an energy diagram shown in Fig. [3j For 
the np-type and n-type interfaces, due to the polar structure of LaA103, the electric field 
through the LaA103 film lifts up the valence band edge of the LaA103 and reduces the 
energy gap. The energy gap of the np-type interface is the energy difference between the 
Ti-d states of the SrTiOs substrate and the O-p states of the SrTi03 capping layer, while 
the energy gap of the n-type interface is given by the energy difference between the Ti-d 
states of the SrTi0 3 substrate and the O-p states on the surface. For the p-type interface, 
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FIG. 3: Schematics of the supercells and energy diagrams for different types of interfaces, a) The 
np-type interface, b) The n-type interface, c) The p-type interface, d) The new quantum well 
systems. 
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since the polarity of LaAlOs is reversed, the electric field through the LaA103 film decreases 
the valence band edge of the LaA103. However, the energy gap of the p-type interface is not 
the most relevant quantity to monitor for the insulating-to-metallic transition. Rather, the 
LaAlOs film reduces the energy difference between the La states on the surface 55] and the 
O p-states of the SrTi03 substrate, which from now on we call the 'La-0 energy difference'. 
In all three cases, the energy gap (np-type and n-type) or La-0 energy difference (p-type) 
diminishes with the LaA103 film thickness, finally disappearing when the insulating-to- 
metallic transition occurs. The minimum number of LaAlOs unit cells necessary to induce 
this phase transition is referred to as the 'critical separation'. 

Our DFT calculations, as well as previous studies [l^Q-lsO], support the above schemat- 
ics. The values of the calculated energy gap versus the LaAlOs thickness are listed in Table 
HT1 for the np-type, n-type and p-type interfaces. The corresponding 'critical separation' 
is the smallest thickness of LaAlOs that makes the interface conducting. The increasing 
thickness of the LaAlOs film reduces the energy gap until it disappears and the interface 
becomes metallic. However, the critical separation depends on the band gap of SrTiOs, 
which is underestimated in DFT calculations (the calculated band gap is 1.85 eV and the 
experimental band gap is 3.2 eV). Taking into account the difference between the calculated 
and experimental band gaps, the realistic critical separation is two more unit cells of LaAlOs 
in addition to the calculated one, resulting in 5, 6 and 8 u.c. for the np-type, n-type and 
p-type interfaces, respectively. 

In order to give a clearer illustration of how the polar field reduces the energy gap, we 
show in the left column of Fig. HI the xy planar average electrostatic potential [40'] and 
the associated macro-averaged smoothed potential along the direction perpendicular to the 
interface for the rip-type geometry. The xy planar average potential is obtained by averaging 
the raw three-dimensional total potential over the interface plane. The associated macro- 
averaged smoothed potential is to convolute the xy planar average potential with a Gaussian 
function using the filting width ~ 0.6A. As the figure shows, the macro-averaged smoothed 
potential in SrTiOs is flat, indicating that there is no internal electric field through the 
SrTiOs (i.e., it is nonpolar). In the LaAlOs layers, the macro-averaged smoothed potential 
is increasing and lifts up the valence band edge of the capping SrTiOs. The energy difference 
between the O-p states of the SrTiOs capping layer and those of the SrTiOs substrate 
increases with increasing LaAlOs film thickness. The right column of Fig. H] shows the 
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FIG. 4: The average potential (left) and density of states (right) of the np- type interface with 
different thicknesses. The upper two panels correspond to the np-type interface with 1 unit cell of 
LaA10"3, and the lower two panels correspond to the np- type interface with 3 unit cells of LaA10"3. 
In the right panels, the red, blue, and orange lines are the total density of states (DOS), the 
atomic projected DOS of Ti-d states, and the atomic projected DOS of O p-states, respectively. 
The vertical green line in the right panels is the Fermi level. 

corresponding density of states (DOS). As the upper panel shows, since the potential increase 
by 1 unit cell of LaAlOs is not large enough to overcome the band gap of SrTiC>3, the system 
is still insulating. In the lower panel, the potential increase by 3 unit cells of LaA103 is larger 
than the band gap of SrTiC>3, so that the Ti-d states in the SrTiC>3 substrate and O-p states 
in the capping SrTiC>3 layer overlap. Then the system becomes metallic. 

Further evidence of the insulating-to-metallic transition is shown by the local density of 
states (LDOS) at the Fermi level of the n-type interface in Fig. \5^l. From the character of 
these states, we can see that we have Ti-d states in the SrTiC>3 substrate and O-p states on 
the surface. Fig. [2] shows the spatial distributions of the conducting electrons and the holes 
integrated over xn-plane at the n-type and p-type interfaces, respectively. The details of 
how to extract out these transferred charge densities are provided in Appendix |A] From Fig. 
[2j it is clear that the conducting electrons and holes behave very differently in the SrTi03 
substrate. The conducting electrons, which occupy Ti-ci states, decay away from the n-type 
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FIG. 5: 3D isovalue surfaces (yellow contour) showing the local density of states at the Fermi 
level for a) the n-type interface and b) the new quantum well systems. The conduction electrons 
(occupying Ti-d orbitals) extend into the SrTiC>3 substrate at the n-type interface, but are localized 
in the embedded T1O2 layer at the QW. Holes occupy the O-p states on the surface in both cases. 

interface over a length scale of ~ 10 unit cells, indicating that the conducting electron gas is 
several nm thick. This result is in qualitative agreement with experimental measurements 
of the superconducting electron gas j^] and other theoretical calculations Q, Q- Unlike 
the conducting electrons, the spatial distribution of the holes resembles that of a particle 
in a box, implying that the holes are very weakly bound to the interface. As elucidated in 
our previous work 281 ]. the reason for this asymmetry is that the conducting electron gas is 
trapped at the n-type interface due to a large interfacial hopping matrix element. We will 
discuss this phenomenon further in Section IIII CI 



2. The quantum well system 

Though the conduction electrons are bound to the n-type interface, they still spread over 
several nanometers. The observed thickness of the electron gas at the n-type interface varies 
widely in experiments, from nanometers to microns, depending on the growth conditions 



44j . Thus it would be desirable to be able to manually control the thickness of the electron 
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gas. Based on the polar catastrophe mechanism, we propose a new class of quantum well 
systems that will allow this functionality. The quantum well system, which is illustrated in 
Fig. EH, has the following composition: SrTi0 3 substrate/2 u.c. LaA10 3 /m u.c. SrTi0 3 
/LaA10 3 capping layer, where m can be varied to control the thickness of the conducting 
electron gas, as we discuss below. Along the (001) direction, this configuration results in 
three interfaces: one p-type interface between the SrTi0 3 substrate and the 2 u.c. LaA10 3 
buffer layer, and two n-type interfaces on both sides of the inserted thin SrTi0 3 film. Both 
of the LaA10 3 thin films are stoichiometric in this sytem, while the embedded SrTi0 3 thin 
film is not. As we have two n-type interfaces, the quantum well systems have one extra Ti02 
layer in the inserted SrTi0 3 film. Like the n-type and no-type interfaces, the two SrTi0 3 
unit cells of the substrate facing the vacuum are fixed to simulate a bulk-like substrate and 
all the other atoms are fully relaxed. 

The reason that this configuration forms a quantum well can also be explained in terms 
of an energy diagram. Fig. [3H shows that the polar fields in the two LaA10 3 thin films point 
in opposite directions. Therefore, the energy of the Ti-<i states in the embedded SrTi0 3 
thin film is lowered relative to the conduction band edge of the SrTi0 3 substrate. (We 
choose two unit cells of LaA10 3 as the buffer layer in order to bring down the Ti-d states 
of the embedded SrTi0 3 film into the band gap of the SrTi0 3 substrate). As the capping 
LaA10 3 layer thickens, the energy of the O-p states eventually becomes higher than that of 
the Ti-d states of the embedded SrTi0 3 film, and hence an insulating-to-metallic transition 
occurs. The difference between the quantum well systems and the n-type interface is that 
once the insulating-to-metallic transition happens, the conduction electrons get trapped in 
the Ti-d states of the embedded SrTi0 3 film rather than those of the SrTi0 3 substrate. 
Therefore, by changing the thickness of the embedded SrTi0 3 film, we can manually control 
the spatial extent of the electron gas while holding the number of transferred electrons 
constant. Consequently, we have two parameters to vary: 1) the thickness of embedded 
SrTi0 3 film to control the spatial extent of electron gas, and 2) the thickness of LaA10 3 
capping layer to control the insulating-to-metallic transition. For simplicity, in the following, 
we focus on one particular subclass in which the embeded SrTi0 3 is narrowed down to a 
single TiC>2 layer. The thickness of the LaA10 3 capping layer will be varied. The quantum 
wells (QW) in the following refer to this particular subclass. 

Like the p-type interface, at this new class of quantum well systems, the energy gap is 
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not the most direct quantity to monitor for the insulating-to-metallic transition. As the 
LaAlOs capping layer is thin, the energy gap of the QW is the energy difference between the 
Ti-d states of the embedded Ti02 layer and the O-p states of the SrTi03 substrate, which 
remains at a constant energy. The more relevant quantity is the energy difference between 
the Ti-d states of the embedded Ti02 layer and the O-p states on the surface. We refer to 
this quantity as the 'Ti-0 energy difference' in the following. Increasing the thickness of the 
LaAlOs capping layer reduces the Ti-0 energy difference until it disappears and electron 
transfer occurs. The minimum thickness of the LaAlOs capping layer that enables this 
electron transfer is defined as the critical separation of the QW systems. 

Table HT1 shows the Ti-0 energy difference and the critical separation of the QW systems. 
The critical separation of the QW coincides with that of the n-type interface, reflecting 
the fact that the band gap of SrTi03 is mainly determined by the Ti-d states and the O-p 
states. Therefore, the absence of Sr atoms does not change the critical separation. The local 
density of states at the Fermi level of the QW are shown in Fig. [5b. Unlike the n-type 
interface, at which the filled Ti ci-states extend over ~10 unit cells in the SrTi03 substrate, 
the only metallic states in the QW systems are localized in the single embedded Ti02 layer, 
as anticipated. 

This new class of quantum wells displays a number of appealing properties that are ab- 
sent at the LaAlOs/SrTiOs n-type interface. First, the thickness of conduction electrons is 
controlled by the inserted SrTiOs film and can be in principle reduced to a single atomic 
Ti02 layer. Second, the electronic properties of the new quantum wells largely depend on 
the polar structure of LaA103. Therefore, the SrTi03 substrate can be replaced by other 
materials on which LaA10 3 can be epitaxially grown, e.g. LSAT (Lao.29Sro.7iAlo.65Tao.35O3). 
This provides more choices of substrates on which to grow this new quantum well structure. 
Finally, the new quantum well serves as a practical way to test the recently proposed hy- 
pothesis that some electrons do not contribute to transport due to Anderson localization. 
Popovic et al. [37J argued that since all the two dimensional states are Anderson localized 
by disorder, the electrons that occupy the lowest Ti-d bands do not conduct and therefore 
the observed sheet carrier density should be much smaller than the 0.5e per two-dimensional 
unit cell predicted by the polar catastrophe mechanism |2^]. Nevertheless, it is not trivial 
to distinguish from which band conduction electrons originate. The thinnest of this new 
type of quantum wells has only one Ti-d band at the Fermi level, corresponding to the sin- 
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FIG. 6: A comparison of the band structures of a) the n-type interface and b) the new QW. 
The n-type interface has 5 unit cells of LaAlOs and the QW has a capping layer of 6 unit cells of 
LaA10"3. The red solid line is the Fermi level. 

gle embedded T1O2 layer, which has a strong two-dimensional Ti-d^ character. Therefore 
it should show strong localization and not contribute to conductivity. By thickening the 
embedded SrTiC>3, one should be changing only the extent of the wave functions and see 
the localization properties. Fig. [6] shows a comparison of the band structures of the n-type 
interface and the QW. These two systems both have approximately O.le per two-dimensional 
unit cell (from Fig. |9]). However, the Fermi level of the n-type interface crosses three Ti-d 
bands, while the QW system has only one Ti-d band at the Fermi level, as expected. 

C. Bound versus unbound carriers at the interfaces 

In this section, we provide a microscopic picture of the quantum states inhabited by 
the carriers at the n-type and p-type interfaces, based on a first-principles extraction of 
tight-binding parameters. As reported in our previous work [28j, the main result is that a 
significant Ti-La interfacial hopping element unique to the n-type interface causes electron 
to bind there. No such mechanism is operative at the p-type interface, explaining the 
derealization of the holes into the SrTiC>3 substrate. As shown below, the Ti-La hopping is 
significantly enhanced at the n-type interface due to the relatively large size of La (i-orbitals 
as well as the spatial proximity of the La and Ti atoms in the neighboring atomic planes. 
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TABLE II: DFT-LDA computed energy gaps and potential changes versus the number of LaA103 
layers i for the np-type, n-type and p-type interfaces (i=0 refers to a pure SrTiC>3 substrate). A 
is the energy gap of the interface systems (np-type and n-type) with the corresponding number 
of LaA103 layers. 5 is the £ La-0 energy difference' (p-type interface) or 'Ti-0 energy difference' 
(quantum well) with the corresponding number of LaAlC>3 layers. V% is the macroscopic potential 
change due to adding the i-th LaA103 layer. Note that 'La-0 energy difference' is not well defined 
when i = (no LaA103 layer). 



np-type 


n-type 


p-type 


QW 


i A(eV) Vi(eV) 


i A(eV) Vi(eV) 


i <5(eV) Vi(eV) 


i J(eV) Vi(eV) 


1.85 


1.85 


N/A 


1.85 


1 1.25 0.60 


1 1.71 0.14 


1 2.59 


1 1.83 0.02 


2 0.55 0.70 


2 1.28 0.43 


2 1.89 0.70 


2 1.30 0.53 


3 metallic 


3 0.57 0.71 


3 1.18 0.71 


3 0.59 0.71 




4 metallic 


4 0.49 0.69 


4 metallic 






5 metallic 





Here we provide an expanded and more detailed explanation and analysis of the hopping 
elements. 

To calculate the on-site and hopping matrix elements, we employ Lowdin atomic or- 
bitals J45] where (r|aR) = (f) a (r — R) is a Lowdin orbital of type a localized around lattice 
position R. However, as the Kohn-Sham wave functions are Bloch states indexed by wave 
vectors k, it is more fruitful to employ Bloch-like superpositions v?a( r ) denned as 



i=yje- k ' R ^(r-R). (2) 



iVk is the total number of ^-points in the Brillouin zone sampling. To extract a tight-binding 
model for the Bloch states at k, we need to calculate the matrix elements ((p^\H\(pp) . Making 
use of the completeness of the Hamiltionian, we obtain: 

(<£\H\<i$ = X>>k)£„ k (nk|^} (3) 

n 

where (r|nk) = ^ n k(r) are the actual Bloch eigenstates of the Hamiltonian H. The overlap 
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(y?k|nk) can be recast easily as follows: 

fa£|nk) = J rfr0;(r)^ nk (r) = (aO|nk) (4) 
Our final operational formula is 

(^|#|^> = ^(a0|nk)£ nk (nk|/30> (5) 

n 

where the overlaps (aO|nk) are automatically computed and reported by the PWscf code. 
Choosing a = (3 gives the on-site energies while a ^ (3 are the hopping elements. The 
orthonormality relation can also be obtained by replacing H with the identity / operator: 

= X>°i nk K nk ^°> = ^ ( 6 ) 

n 

Eel-© is a good criterion to check the truncation in the infinite summation over the band 
index n. In our calculations, we include bands with energies up to 29 eV above the Fermi 
level, so that K^IVa)! > 0-99 f° r au Lowdin orbitals considered and |(<^a|y^)| < 5 x 1CT 4 
for all the pairs of two different atomic orbitals. 

At the n-type interface, direct projection of the Bloch states onto the atomic orbitals 
shows that the character of the lowest band accommodating the conduction electrons is 
mainly Ti-d xy with a small component of La.-d xy . The minimum of these occupied bands 
is at r (k = 0). Therefore we build our tight-binding model on the subspace composed of 
Ti-d xy and La-d xy orbitals and calculate the following on-site and hopping matrix elements: 



#oo = ( ( Ph a ,-d my \ H \ c Ph a ,-d xy ) ( 7 ) 

H J3 = (^-djHl^-dJ ( 8 ) 

Ho, = H; o =(^_ d jH\^ TiJ _ d J (9) 

Hij = H*, = & T Tii _ dxy \H\^_ d J (10) 



where 'La' in Eq. (j7j) and Eq. ([9]) is the La atom in the LaA103 layer at the n-type interface. 
i (or j) labels the Ti atoms in the SrTiC>3 where i — 1 (or j — 1) is in the very Ti02 layer at 
the rz-type interface, and increasing i (or j) refers to the Ti02 layers that are further away 
from the interface. 
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FIG. 7: The on-site matrix elements of the n-type interface. The squares are the on-site matrix 
elements. Site is the La atom. Sites 1-6 are the Ti atoms. The triangles are the 'new' onsite 
matrix elements after taking into account the hopping effects. The data can be read from Table 
IIII1 a is the lattice constant of SrTiC>3. a) The n-type interface with 2 unit cells of LaAlC>3. b) 
The n-type interface with 4 unit cells of LaA103. The purple arrows illustrate the effect of Ti-La 
hopping as per Eq. (fT2|) and Eq. (fl~3j) . 



The n-type interface with 2 u.c. of LaAlOa is insulating (before the polar catastrophe) 
and that with 4 u.c. of LaA10 3 is conducting (after the polar catastrophe). Table Hill shows 
the on-site and hopping matrix elements of the n-type interface before and after the polar 
catastrophe, respectively. The on-site matrix elements are also shown in Fig. [7J We can see 
that the polar catastrophe does not result in significant difference in either the on-site or the 
hopping matrix elements. Therefore, in the following we mainly focus on the tight-binding 
Hamiltonian after the polar catastrophe, which is more relevant to describe the conduction 
electron gas. 

Fig. [7J and Table II I II show that the on-site matrix element of the first Ti atom (closest to 
the interface) is not the lowest one before or after the polar catastrophe; instead, the second 
Ti atom has the lowest potential. Therefore, one would expect to find the highest electron 
density on the second Ti atom, in direct contrast to what is actually found in Fig. |2fc. 
However, there is a very large hopping matrix element between the La and Ti atoms at the 
interface, at least 100 times larger than all the other hopping matrix elements between Ti 
atoms. This is due to two factors: one is that La-5<i orbitals are more extended than Ti-3<i 
orbitals, which leads to a larger overlap, and the other is that the distance between La and Ti 
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atom at the n-type interface is a/3/2 times the distance between neighboring Ti atoms 
Consequently, to the leading order [56|, we neglect the hopping between adjacent Ti atoms 
and other higher order hopping matrix elements. The Hamiltonian becomes block- diagonal 
and the only block we need to diagonalize is the leading 2x2 Ti-La submatrix right at the 
interface: 

h = ( H °° H A (ii) 

\Hio Hu J 

Diagonalization of Eq. f lTT]) with the tabulated values gives the following eigenvalues and 
eigenvectors: 

e x = 12.99eV, |6) = 0.98|La-4 y ) - O.lSlTi 1 -^) (12) 

e 2 = 9.04eV, |f 2 ) = 0.18|La-4 y ) + O^lTi 1 -^) (13) 

After diagonalization, we find a pair of bonding and anti-bonding states at the n-type 
interface. The energy of the binding state e 2 is even lower than H22, showing that the 
electrons prefer the Ti site right at the interface. Therefore, after the insulating-to-metallic 
transition occurs, most electrons reside in this bonding state, which leads to an electron gas 
bound to the interface. In other words, the large Ti-La hopping matrix element is critical in 
creating a strongly bound state for the electrons at the interface. Moreover, this hopping is 
intrinsic to the n-type interface geometry as it does not significantly change its magnitude 
before and after the polar catastrophe. This also explains the overall similarity of the electron 
distributions at the symmetric double n-type superlattice (Fig. [2^) and the stoichiometric 
n-type interface (Fig. [2b). Despite differing in the lack or presence of oppositely charged 
carriers and polar fields, the two systems share the same Ti-La hopping and thus the same 
binding force to the interface. 

On the other hand, at the p-type interface, holes are found to occupy bands with essen- 
tially pure O-p characters and the maximum of those bands is located at M (k = (-, 0)). 
The tight-binding model is based on the following matrix elements: 

( 14 ) 

where i,j refer to different O atoms and /1, v = x,y,z refer to different O-p orbitals. We 
calculate all these matrix elements at the p-type interface before and after the polar catas- 
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trophe. We find that the hopping elements between nearest neighbor O atoms are all close 
to 0.6 eV, that there is no order of magnitude difference in the hopping elements, and that 
there is no discontinuity in going across the p-type interface. This is no surprise because 
the O atoms form a continuous network across the p-type interface and the hopping matrix 
elements do not sensitively depend on the nature of the surrounding cations. Therefore, the 
only difference we expect from the symmetric double p-type results (compare Fig. [2b and 
Fig. [2H) is that the presence of the polar field at the p-type interface will drive the holes 
from the LaAlOs into the SrTiC>3. However, once in the SrTiC>3, the holes do not feel any 
strong preference for the interface and therefore diffuse further into the SrTiC>3 substrate. 
As our calculations have finite film thicknesses, we expect that the profile of holes resembles 
the lowest state of a free particle in a box and thus have a maximum density in the middle 
of the film, as borne out by Fig. [2U. 



D. Oxygen vacancies repulsion from the interfaces 

The theoretical results we have discussed to this point have concerned ideal interfaces 
with sharp boundaries and lxl in-plane periodicity. An important question is whether 
these assumptions are reasonable for describing the experimentally realized interfaces. This 
question is obviously very difficult to answer in general as the number of possible ways of 
perturbing the ideal interfaces is enormous and includes intermixtures, impurities, vacancies, 
interstitials, off-stoichiometry, etc. Below, we focus on answering one simple but important 
question for the most prevalent type of imperfections: do oxygen vacancies present in the 
SrTiC>3 substrate have a strong preference for the SrTiC>3 region close to the interface? Based 
on both our first principles results and general physical considerations, the answer appears 
to be no. As we show below, the interface SrTiC>3 region should be relatively free of oxygen 
vacancies compared to the bulk SrTiC>3 substrate. 

Our approach is motivated by experiments that infer a concentration of one oxygen 



vacancy per four two-dimensional unit cells at the p-type interface [27|]. If one assumes 
each oxygen vacancy donates fully two electrons, this is precisely the amount required to 
compensate the 0.5 holes per two-dimensional unit cell from the polar catastrophe. The 
same experiments are also used to infer a non-zero (but smaller) concentration of oxygen 
vacancies at the n-type interface j^. To determine the most energetically favorable position 
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TABLE III: On-site and hopping matrix elements of the n-type interface before and after the polar 
catastrophe, respectively. I is the thickness of LaA103. I = 2 and 4 unit cells (u.c.) correspond 
to before and after the polar catastrophe, respectively. We list all the nearest neighbor hopping 
matrix elements and the largest next nearest neighbor hopping matrix element. The first and third 
columns are the on-site matrix elements of Ti d-states and hopping matrix elements before the 
polar catastrophe. The second and fourth columns are the on-site matrix elements of Ti d-states 
and hopping matrix elements after the polar catastrophe. 



H u (eV) 


H l3 (meV) 


i 


I = 2 u.c. 


I = 4 u.c. 


ij 


I = 2 u.c. 


I = 4 u.c. 





12.99 


12.87 


1 


-697 


-696 


1 


9.13 


9.18 


1 2 


-6.39 


-1.3 


2 


9.09 


9.14 


2 3 


-6.80 


-0.2 


3 


9.10 


9.15 


3 4 


-7.03 


-0.3 


4 


9.10 


9.19 


4 5 


-7.07 


-0.3 


5 


9.10 


9.16 


5 6 


-7.07 


-0.4 


6 


9.10 


9.17 


6 7 


-7.05 


-0.4 


7 


9.10 


9.18 


7 8 


-7.03 


-0.5 


8 


9.10 


9.19 


8 9 


-6.98 


-0.5 


9 


9.10 


9.19 


2 


-0.87 


-7.2 



for oxygen vacancies, we compute the formation energy of one monolayer of oxygen vacancies 
at the same areal density (1/4 per two-dimensional unit cell) at various positions inside the 
SrTiC>3 close to the interface with LaAlCV We perform calculations on both p-type and 
n-type interfaces using the unit cells shown in Fig. |HJ 

Our simulation cells include a 2 x 2 x 5 SrTiC>3 film and a 2 x 2 x 1 LaA103 film and ~ 
25 A of vacuum. The substrate termination is SrO in order to minimize any surface effects 



281 ] . The SrTiC>3 unit cell next to the vacuum is fixed to simulate the bulk-like substrate 



and all the other atoms are relaxed. For a reference, we compute the formation energy 
of an isolated oxygen vacancy in bulk SrTiC>3 in a 2 x 2 x 5 supercell with one oxygen 
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vacancy. Starting from the interface, we place the oxygen vacancy in the SrO or Ti02 layer 
of the first three SrTiO"3 unit cells and compute the formation energy (see Fig. [8] for the 
numbering nomenclature of the layers). We show the results in Table HVl As the oxygen 
vacancies move away from the p-type interface, the formation energy decreases, approaching 
a constant value of 5.4-5.5 eV in the bulk-like regions of the SrTi03 films. The computed 
bulk formation energy is slightly lower (5.18 eV); presumably the small difference is due to 
the fact that the SrTiO"3 films in the interface geometry are rather thin (only 5 unit cells) 
and not yet fully in the substrate limit. Nevertheless the energetic trend clearly shows the 
interface region repels the vacancies. 

Our results show that there is no energetically favorable binding of oxygen vacancies to 
either the p-type or ra-type interface. A simple physical picture can explain the repulsion 
of the oxygen vacancies from both types of interface. An isolated, neutral oxygen vacancy 
in SrTiC>3 is a donor that binds two electrons (O has formal charge 2 ~). The electrons 
reside on localized states composed of the <i-orbitals of the vacancy's surrounding Ti atoms; 
the energy of these states is close to the SrTiC>3 conduction band edge. Thus, oxygen va- 
cancies can be described as a type of hydrogenic system with bound electrons. As the 
vacancy approaches the interface, the bound electrons experience what are essentially hard 
wall boundary conditions since the conduction band edge of LaA103 is ~ 2 eV higher than 
that of SrTiCV As the vacancy approaches the hard wall, its energy increases due to the 
electron confinement effect [46| . Given the basic physical principles behind this argument, it 
is clear that it applies to both interfaces. Futhermore, it suggests that even though our first 
principles results are for relatively high densities of oxygen vacancies in a particular con- 
figuration, lower densities and more positionally disordered oxygen vacancies in the SrTiO"3 
will still be repelled from the interface region. 

The fact that oxygen vacancies are repelled from the interfaces provides a self-consistent 
picture for our computations of ideal interfaces. Namely, the dominant defects (oxygen 
vacancies) should not be present in the immediate vicinity of the interfaces. For example, 
for the p-type interface, our theory predicts extended band states for the holes that diffuse 
substantially into the SrTiC>3 substrate. Since the vacancies are themselves in the bulk of the 
substrate, they can trap the holes, rendering them immobile. For a relatively thick SrTiO"3 
substrate, this situation is likely since only one oxygen vacancy per four two-dimensional 
unit cells, distributed over the three-dimensional volume of the substrate, is required to 
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provide the compensating number of electrons. 

We would like to comment that experimentally Nakagawa et al. j^] found 32 ± 6% of 
oxygen vacancies per two-dimensional unit cell at the p-type interface. Their conclusion is 
based on a least-square fit of EELS spectra to the reference spectra of bulk SrTi03, bulk 
LaAlOa and oxygen-deficient SrTi03_,5 with 6 = 1/4. The discrepancy between our results 
and those experiments may be due to the following two non-mutually exclusive reasons. One 
is that in our simulations, we include only one unit cell of LaA103 owing to the computational 
limitations. More realistic simulations would include thicker LaA103 films to simulate the 
system after the polar catastrophe when carriers appear in the vicinity of the interfaces. 
The presence of carriers might in principle affect the formation energy of oxygen vacancies, 
which needs to be checked in much larger calculations in future work. The other reason 
is the interpretation of EELS spectra in experiment. The experimental fitting is based on 
the assumption that any deviation of O-K edge EELS from the bulk SrTiC>3 is caused by 
the presence of oxygen vacancies. However, even without oxygen vacancies at the p-type 
interface, we observe in our calculations that the atoms close to the interface move away 
from their bulk positions due to the broken symmetry at the interface. The EELS spectra of 
the distorted SrTiC>3 could be different from that of the reference bulk structure. Therefore, 
a more accurate (and more difficult) fit would take into account the deviation of the EELS 
spectra of distorted SrTiC>3 from the bulk and other possible factors which could also affect 
the EELS spectra, such as intermixture and nonstoichiometry. 



E. Thickness dependence of sheet carrier density 

In experiment, the sheet carrier density was observed to display a thickness dependence, 
but a range of results exist. The sheet carrier density depends on both sample growth 



conditions and post-annealing [441 ] . In this section, we present first-principles calculations of 
the sheet carrier density for different LaAlC>3 film thicknesses. We study the np-type, n-type 
and QW systems, and discuss the similarities between these different interface systems. 

We define the 'sheet carrier density' as the sum of all the conduction electrons per square 
unit cell in the Ti ci-states of the SrTiC>3 substrate (for n-type and np-type interfaces) or 
in the Ti d-states of the single embedded TiC>2 layer (for the QW systems). Fig. |9] shows 
the sheet carrier density at different thicknesses for all three types of interface. The results 
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FIG. 8: Schematics of the simulation cell for 4 monolayer oxygen vacancy calculations. One oxygen 
vacancy is placed in either SrO and TiO"2 of the first, second and third layer, respectively, a) The 
p-type interface, b) The n-type interface. 

of n-type interface are similar to what Son et al. recently obtained [30 1 . Below the critical 
separation in each system, there are no mobile carriers since the system is insulating. Above 
the critical separation, the sheet carrier densities increase in all three cases, following a 
similar manner. In order to further reveal the similarities among these interface systems, we 
show in Fig. [10] the relation of the sheet carrier density a s versus the internal electric field 
E through the LaAlOs. The internal electric field is determined as follows. We calculate the 
macro-averaged smoothed total potential along the (001) direction. The internal electric field 
is the slope of the smoothed potential. Though the geometry details differ, the figure shows 
that the (E, a s ) relations are almost identical for all three interface types. This suggests 
that a single continuous model will be able to describe the behavior of all three interfaces. 
We discuss such a model in Section IIV A[ 

We note that in our calculations of the sheet carrier density, we include all of the electrons 
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TABLE IV: Formation energies E of \ monolayer of oxygen vacancy at the p-type and n-type 
interface, j monolayer of oxygen vacancy is in either zth SrO atomic layer or ith Ti02 atomic 
layer. For reference, the formation energy in the bulk is 5.18 eV. 



p-type 


n-type 


position of oxygen vacancy 


E (eV) 


position of oxygen vacancy 


E (eV) 


1st SrO layer 


6.25 


1st Ti02 layer 


5.65 


1st TiC>2 layer 


5.85 


1st SrO layer 


5.65 


2nd SrO layer 


5.69 


2nd Ti02 layer 


5.56 


2nd Ti02 layer 


5.56 


2nd SrO layer 


5.46 


3rd SrO layer 


5.60 


3rd Ti02 layer 


5.57 


3rd Ti02 layer 


5.50 


3rd SrO layer 


5.41 



in the conduction bands, regardless of whether they in fact contribute to the observed 
conductivity or not. We find that after the critical separation, the sheet carrier density 



film, asymptoting to 0.5 e/a 2 . This 
23]. but conflicts with others jj], in 



increases continuously with the thickness of the LaAlOs 
behavior qualitatively agrees with some experiments p, 
which the sheet carrier density remains almost constant (15 times smaller than 0.5 e/a 2 ) 
after the critical separation. It is possible that the low sheet carrier density in Ref. j4j may 



be explained in terms of Anderson localization induced by disorder [30|, |37|. 



F. External Field Effect 



We have shown that there is an insulating-to-metallic transition as a function of LaA103 
thickness both in the SrTiOs/LaA103 interface and the quantum well (QW) systems. In 
this section, we examine the effects of inducing this transition at subcritical thicknesses via 
an applied electric field. We focus on the n-type interface and the QW system with 3 u.c. 
of LaA103 (for both systems, the critical separation is 4 u.c. in the DFT simulations). To 
study the field effect, we apply a homogeneous electric field perpendicular to the interface. 
In order to avoid artificial effects from periodic boundary conditions (see Appendix IBI) . 
we apply two electric fields of equal magnitude and opposite direction to each half of a 
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FIG. 9: Comparison of the sheet carrier density computed with the continuous model (see Section 
IIV Ap and DFT simulations. The red, orange and green lines are the model calculations of the 
QW, the n-type and the np-type interfaces, respectively. The squares, triangles and circles are the 
DFT results of the QW, the n-type and the np-type interfaces, respectively. 

symmetric simulation cell (see Fig. HTjl . Since the system is mirror symmetric, for simplicity 
we only focus on half of the simulation cell. At the n-type interface, the external electric field 
is parallel to the intrinsic electric field through the LaA103 film, and in the QW systems, 
the direction is parallel to the field through the LaA103 capping layer. Based on the polar 
catastrophe picture, such external electric fields will add to intrinsic electric fields and further 
reduce the energy gap, so that the insulating-to-metallic transition will occur at a subcritical 
separation. 

We performed the simulations with different magnitudes of external electric field on both 
the n-type interface and QW systems. We need to elucidate a subtle point here: in the QW, 
the external electric field reduces the energy difference between the Ti d-states of the embed- 
ded Ti02 layer and O p-states on the surface. This is the 'Ti-0 energy difference' introduced 
in Section IHIB 21 Once this energy difference disappears, the insulating-to-metallic transi- 
tion occurs. This energy difference is well denned because we can identify two Bloch states, 
one which characterizes the Ti ci-states of the embedded TiC>2 layer and the other which 
characterizes the O p-states on the surface, and simply take their energy difference. At the 
n-type interface, the external electric field also reduces the energy difference between the Ti 
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FIG. 10: The sheet carrier density versus the internal electric field through the LaA10"3 thin film. 
The red squares, orange triangles, and green circles are for the quantum well, the ra-type, and the 
rap-type interfaces, rspectively. The blue line is the result of the continuous model described in 
Section ESI- 

d-states of the first Ti atom at the interface and the O p-states on the surface. However, 
this energy difference is not easily identified because the external electric field also exists 
throughout the SrTiC>3 substrate, so the Bloch states that accept the transferred electrons 
are now composed of a mixture of all the Ti d-states in the SrTiC>3 thin film. For the 
purposes of this section, we simply extract the energy gap of this ra-type interface system 
(between the O p-states on the surface and the lowest occupied Ti d-states) versus the ex- 
ternal field. In Section IIVBI we will discuss more fully the relation of this energy gap and 
what would happen in an actual experiment. 

In Fig. [12] we show the computed energy gap of the ra-type interface and the Ti-0 energy 
difference of the QW as a function of external electric field. As the figure shows, both 
monotonically decrease with increasing external electric field. Thus, the DFT simulations 
demonstrate how an external electric field can induce an insulating-to-metal transition. 

IV. CONTINUOUS MODEL 

In order to shed more light on the nature of the LaAlOa/SrTiOs interface, we develop a 
simple model which approximates the LaAlOs film as a homogeneous continuous medium. 
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QW 




FIG. 11: The supercells used for the field effect calculations. For the QW and the n-type interface, 
both simulation cells are mirror-symmetric. The external electric field is also mirror-symmetric. 
The dashed line in the middle is the symmetry axis. The embedded TiO"2 layer in the QW and the 
n-type interface LaO/TiO"2 are also highlighted by the dashed line. 



FIG. 12: Energy differences versus the external electric field. The red squares and orange triangles 
correspond to the DFT-computed Ti-0 energy difference for the QW systems and the energy gap of 
the n-type interface, respectively. The solid lines are the results of the continuous model described 
in Section IIVBI 

This model shows that the polar and dielectric properties of LaAlOs largely determine the 
thickness dependence of the sheet carrier density and the external field effect behavior. 
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FIG. 13: Schematic illustration of the continuous model. In the ionic limit, the LaAlC>3 thin film 
can be considered as a serial connection of capacitors. The electric field through LaA103 is like an 
impulse and the resulting potential takes a stair shape. The dashed purple lines are the average 
electric field and the average potential. 

A. Thickness dependence of sheet carrier density 

In this section, we use model calculations, aided by DFT results, to quantitatively describe 
the insulating-to-metallic transition that occurs at the n-type and np-type interfaces and the 
quantum well (QW) systems. As detailed below, we find that the electric field dependence 
of €l is critical in obtaining quantitative accuracy. 

Before the insulating-to-metallic transition occurs, the alternating positively charged 
(LaO) + and negatively charged (AIO2) - layers in the LaAlOs film can be idealized as a 
series of capacitors (see Fig. [131). As the figure illustrates, the thickness of each capacitor is 
only half of the unit cell of LaAlOs (there is zero voltage drop across the other half). There- 
fore the inner electric field through the LaAlOs behaves like an impulse and the resulting 
potential takes the shape of a stair function. If we approximate the LaA103 thin film as a 
continuous homogeneous media and average this electric field, we can determine the average 
internal electric field: 

E=l^ = ^ } (15) 
2 ei e L 

where a = e/a 2 is an electron per two-dimensional unit cell, is the dielectric constant of 
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LaAlC>3, and 



a = Y = 2^- (16) 

Once the insulating-to-metallic transition occurs, some electrons transfer across the interface 
to counteract the internal electric field through the LaA103 thin film. The new average 
internal electric field through the LaA103 film is then given by: 

E= A ^-^\ (17) 

where a s is the transferred electron density, or sheet carrier density. It is easy to recast 
Eq. ([I7D into: 

a s (E) = a-^ (18) 

where a = 0.5e/a 2 . We also take into account that the dielectric constant depends on the 
internal electric field E. The field dependence can be phenomenologically described by the 



Landau theory, and takes the following approximate form [43, |48j (also see Appendix ICl) : 




(19) 

where eo is the dielectric constant of LaAlOs at vanishing electric field and So is a charac- 
teristic electric field. Eq. ffl8|) and Eq. ffl9|) establish a relation between the internal electric 
field E through the LaA103 and the sheet carrier density a s , both of which can be calculated 
independently from DFT simulations. We also performed a separate slab calculation (see 
Appendix IC II for details) and determined the parameters eo and So in Eq. ( fT9l) as: 

e = 40.95, £ = 0.15V/A (20) 

Fig. [TUl compares (E,a s ) computed via DFT and the model. We can see from Fig. [TOlthat 
(E, a s ) almost lies on the same curve in all three cases, which does not depend on the details 
of the structure. The continuous model shows a good agreement with the DFT results, 
and as the internal electric field decreases (i.e. the thickness of the LaA103 increases), the 
continuous model becomes almost exact. This is expected because in the limit of an infinitely 
thick LaA103 film, the interface and surface regions become neglegible. Our results show 
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that even with only 7 u.c. of LaA103 (i.e. the internal electric field < 0.12 V/A), the 
continuous model already works very well. 

However, it is not easy to directly measure the internal electric field through the LaA103 
experimentally because the LaAlOs film is only a few unit cells thick. It is more useful to 
relate the sheet carrier density to the nominal thickness of the LaAlOs film (i.e. the number 
of LaA103 unit cells grown on top of the SrTi03 substrate). Therefore, we need to find how 
the internal electric field through the LaAlOs depends on the LaAK^ thickness. 

When the polar catastrophe takes place at the n-type interface, electrons leave the surface 
O-p states and are transfered into the Ti-d states. The charge transfer is halted once a 
common Fermi level is reached. Approximately, this means that the Ti-<i conduction band 
edge of the SrTiOs substrate is at the same energy as the O-p surface valence band edge. 
Looking back at Fig. [31 we can see that there is a remanent field through the LaAlOs film 
and thus a potential difference. The existence of the electric field is also visible in Fig. [TOl 
from our calculations. The potential difference across the LaAlOs thin film, denoted by K, 
can be related to the remanent internal electric field by: 



where d is the thickness of LaAlOs. In principle K could depend on the thickness of LaAlOs. 
But the detailed calculations show (see Fig. IT^l) that K quickly saturates to a constant as the 
thickness d gets larger. This could be understood because K is essentially determined by the 
energies of the relevant electronic states and the band offsets, which are interface properties 
and have little dependence on the thickness of LaAlOs. We therefore approximate K as 
a constant. We illustrate later that this approximation does not significantly change the 
physics in the continuous model. 

Inserting Eq. (12T1) into Eq. (f!8l) and Eq. ( 1191) . we eliminate E and obtain: 



Eq. f l2"2"|) highlights some qualitative features of the sheet carrier density. The critical sepa- 
ration is the smallest d which makes a s > 0. Below the critical separation, the sheet carrier 
density is zero. Above the critical separation, the sheet carrier density gradually increases 
and then saturates to a = 0.5e/a 2 in the limit of large d. 



eEd = K 



(21) 




(22) 
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FIG. 14: The potential difference across the LaA10"3 thin film, defined by eEd, of the quantum 
well (QW), the n-type and the np-type interfaces. Squares, triangles, and circles correspond to the 
QW, the n-type, and the np-type interfaces, respectively. E is the average internal electric field 
along the (001) direction through the LaAlOs thin film and d is the nominal number of LaA10"3 
unit cells (note that the SrTi03-strained LaAlG^ lattice constant olao = 0.953osto)- The dashed 
lines highlight the data point of the largest thickness in each interface configuration. 

The parameter K in Eq. ( 122]) is determined as follows. We simulate n-type, np-type 
interfaces and QW with different thicknesses and calculate the macro-averaged smoothed 
electric field through the LaA103 at each thickness. In Fig. [TU we show the potential 
difference across the LaAlC^ thin film, defined by eEd. We can see that in all three cases, eEd 
quickly saturates as d gets larger. Since we do not have a microscopic model to describe how 
the potential difference eEd changes with the thickness d, we use the last data (corresponding 
to the largest d) to determine the values of K, which are highlighted by the dashed line in 
Fig. [TU In each case, we find that 



QW 



3.16 eV, K„ 



2.64 eV, K np 



2.40 eV. 



(23) 



The K in all three cases turn out to be different from each other and larger than the DFT- 
computed band gap of SrTiC>3 (1.85 eV). This is partly because the band offsets of these three 



13 



15j can not explain 



cases are not the same. However, the small valence band offsets alone 
the big difference between the K value and SrTiC>3 band gap. A more important factor is 
that we treat E as a uniform electric field throughout the LaAlOs thin film. However, close 
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to the interface and surface, the internal electric field should be different from that in the 
middle of the LaAlOs film but difficult to average in those regions. Therefore we use an 
approximate relation (K = eEd) to determine the potential difference, but as shown later, 
this approximation is good enough to reproduce the DFT-calculated thickness dependence 
of sheet carrier density. 

Equipped with K, e and £ , we can calculate the sheet carrier density at different 
thicknesses using Eq. f l22|) . We compare the sheet carrier densities calculated from the 
continuous model with those from the DFT simulations. The results are shown in Fig. |9j 
The agreement is good in all three cases even though Eq. (f2Tl) is not exact. 

We would like to comment that recently Son et al. [30] also calculated the thickness 
dependence of sheet carrier density and they made a similar model by assuming that the 
dielectric constant of LaA103 has no field dependence and is simply constant. Rather than 
Eq.( )22|) . they had a simpler formula 



where a and d c are two constants determined by fitting. They found a = 0.455e/a 2 , close 
to the ideal value 0.5e/a 2 . Our analysis is different in that we fix a = 0.5e/a 2 , which is 
necessitated by the polar catastrophe mechanism, but take into account that the dielectric 
constant of LaA103 has a strong field dependence, as confirmed by the DFT calculations. 

B. External electric field 

In order to get a more quantitative description of the external field effect, we apply 
the continuous model established in the previous section. In principle, at the subcritical 
separation (3 u.c), the continuous model should break down because the LaAlOs film is so 
thin that the interface and surface regions cannot be neglected and the LaAlOs is no longer 
uniform. However, we simplify the situation by assuming that LaA103 is still a homogeneous 
media with the dielectric constant given by Eq. f|T9l with the parameters Eq. fl20l . but has 
an effective thickness d e ff. In order to simplify our notations, we use E instead of E to 
denote the averaged macroscropic fields. 

Based on the above argument, for the quantum well systems we have: 




(24) 
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A = A - (E L - E%)d L eff , (25) 

where l U stands for LaA103. A is the Ti-0 energy difference at a given E ext and E L is the 
internal electric field. Ao and Eq indicate the values of A and E L , respectively, at vanishing 
external electric field: 

E ~^wv (26) 

and 

Equations (I25 [ |2"6" | and l2"Tj) establish the relation between A and E ext . The effective thickness 
cannot be given a priori from the model. Instead, we fit the data and find d\^ = 5.9A = 1.6 
u.c. A comparison of (E ext , A) for the quantum well systems using DFT calculations and 
the model calculations is shown in Fig. [12j 

For the n-type interface, we also need to account for the potential drop in the SrTiC>3 
due to the external electric field. Therefore we have: 



A = A - (E L - E^)d L eff - E b d b eff (28) 

where £ 5" stands for SrTiC>3. A is the energy gap of the system and Ao is the value of A at 
vanishing external electric field. Note that in the absence of an external electric field, the 
SrTiC>3 film is unpolarized in the continuous model. E s and E ext are related by: 

E>=*^ (29) 



We use the Berry phase method [49(] to calculate es{E) and fit the data using Eq. (Tl9|) (see 



Appendix IC 21 for details). The results of the fit are: 

4 = 309.6, 8% = 49.2V//im = 4.92 x 1(T 3 V/A (30) 

We find that the values of efe^ in different systems are different. Hence, in principle we have 
two fitting parameters in Eq. fl28l : d^j and c?f//- However, it turns out that Eq. fl28|) very 
insensitively depends on d s e ^. This is expected because the dielectric constant of SrTi0 3 is 
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at least 10 times larger than that of LaA103 and the nominal thicknesses of LaA10 3 and 
SrTiOs are both 3 u.c. in the DFT calculations (see Fig. ITT]) , and thus the potential drop 
across the SrTi03 is much smaller than that across the LaA103 film. Therefore we just use 
the real thickness of SrTiOs d^jff = 11.1 A (from the relaxed structure) and fit the data with 
d\ n . We find: d% ff = 7.5 A~ 2 u.c. Fig. QJ compares the values of (E ext , A) of the n-type 
interface computed by DFT and the model, respectively. 

Fig. [J_2] shows that the critical electric fields E c ext extracted from the DFT calculations 
for the n-type interface and the QW are 0.79 and 1.15 V/A, respectively. In the continuous 
model, the critical voltage across the whole sample is given by: 



V B = d LAO + -^sto (31) 



Et 

"LAO H 

CLAO CSTO 

where g?lao and g?sto are the thickness of LaAlOs and SrTiOs, respectively, a = 0.5e/a 2 
and E c ext is the critical external electric field. In order to estimate the critical voltage for the 
experimental setup, we need to use the thickness of the SrTiOs substrate in experiment. 

Here we need to clarify some subtleties: in the presence of applied external field, DFT 
simulations do not realistically represent the spatial distribution of conduction electrons. 
Fig. US] illustrates the reason for the discrepancy. In the DFT simulations, the conduction 
electrons fill the lowest energy states. Since the "external" field bends the SrTiOs conduction 
bands, these states are located at the bottom surface of the SrTi0 3 substrate, as shown 
schematically in Fig. [15a . Although in principle the experimental situation is the same, 
realistically, electrons can only tunnel a few unit cells from the surface; thus they get trapped 
at the n-type interface, as illustrated in Fig. [ToT d. The trapping mechanism is the Ti-La 
interfacial hopping, which lowers the energy of the Ti atom closest to the interface and 
creates the energy barrier through which the electrons must tunnel. Hence, the real critical 
field has to be larger than the computed one, so that the energy of surface states equals 
that of the interface states (Fig. [T5b). The critical external field E c ext we obtain from DFT 
simulations is therefore only a lower bound. 

Now we make a simple estimation of the lower bound: at low temperatures (< 5 K), 



the dielectric constant of SrTi0 3 can be as large as 2.5 x 10 4 50| and the thickness of the 
SrTi0 3 substrate is typically ~ 1 mm. Considering that 6lao — 30 and g?lao — 10 A, the 
second term in Eq. (l3TT) is dominant. Therefore, for the n-type interface, we have: 
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FIG. 15: a) The spatial distribution of conduction electrons in the DFT simulations. The conduc- 
tion electrons occupy the lowest energy states available, which are located at the bottom surface, 
b) The spatial distribution of conduction electrons in the actual experiment. The conduction 
electrons get trapped at the n-type interface by the tunneling barrier due to the self-consistent 
potential well formed at the interface that largely stems from the Ti-La hopping. The length of 
arrows illustrates the magnitude of external field. 



yn-typc _ _k rfgT0 ~ 32 V (32) 

csto 

and for the QW, 

V C QW ~ — d STO ^ 480 V (33) 

^STO 

As these estimates are based on data from DFT simulations in which the band gap of SrTiC>3 
is underestimated, they provide only lower bounds on the critical voltage. Experimentally, 
the critical voltage of the n-type interface with a 1 mm thick SrTiC>3 substrate has been 
measured to be ~ 70 V at low temperature 4j, which is much smaller than our estimates. The 
discrepancy may imply that the simple model of an ideally sharp interface is not sufficient to 
model the external field doping, due to additional details that play an important role at the 
n-type interface. Various types of atomic reconstructions, such as cation intermixtures at 



the interface 



27j or oxygen vacancies on the surface may possibly trigger an insulating- 



to-metallic transition, leading to a smaller critical voltage. 
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V. CONCLUSION 



We present a detailed study of ideal LaAlOs/SrTiOs interfaces and a new class of quan- 
tum wells. In both systems, we confirm an intrinsic insulating-to-metallic transition. The 
observed transition can be triggered either by thickening the LaAlC>3 layers or by imposing 
an external electric field, both of which can be explained by the polar catastrophe mecha- 
nism. We show that for both the n-type interface and QW, the realistic critical separation 
(taking into account the underestimation of DFT band gap) is around 6 unit cells. We 
also show that, given a typical SrTiC>3 substrate thickness of 1mm, the lower bound for the 
critical voltage necessary to induce an insulating-to-metallic transition is estimated to be 
~ 300 V for the ideal n-type interface and ~ 500 V for the QW. In addition to theoretically 
demonstrating the observed physical properties, we provide a microscropic explanation of 
the observed binding of conduction electrons at the n-type interface, a phenomenon which 
can not be described by the polar catastrophe mechanism alone. We demonstrate that the 
large La-Ti hopping matrix element at the n-type interface, which is absent in both bulk 
constituents and is unique to the n-type interface, lowers the energy of the Ti atom at the 
interface relative to all other Ti sites, thus binding the electrons. Futhermore, we develop a 
continuous model that captures the essence of polar structure and dielectric properties and 
predicts the thickness dependence of sheet carrier densities. 

We conclude by discussing some of the outstanding issues related to this system. Our 
understanding of the atomically sharp interfaces and QW suggests that the electronic prop- 
erties of LaAlOs/SrTiOs interfaces observed in experiments are unlikely to be solely a result 
of electronic reconstructions based on the polar catastrophe mechanism. For example, the 
theoretical critical separation, predicted by the polar catastrophe picture, is 6 u.c 
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301 ]. which is larger than the experimental value (4 u.c.) [4j. In addition, while the calculated 



thickness dependence of the sheet carrier density agrees qualitatively with some experiments 



15, 



l23j ] in that the sheet carrier density increases with LaAlOs thickness, it disagrees with 
other experiments {4 ] which show an almost constant sheet carrier density above the critical 
separation. The sheet carrier density we calculate here includes all the conduction electrons. 
Whether some of them may be prone to the Anderson localization [37] is an issue beyond 
;he scope of the DFT calculations. However, the possibility of multiple type of carriers 



5l| is an interesting question that warrants further study in both theory and experiment. 
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Another remaining puzzle is that the polar catastrophe predicts the coexistence of electrons 
at the interface and holes one the surface, while experiments find that the surface region 
is insulating and only electron-like carriers are detected in the transport measurement ^]. 
Recently, x-ray photoemssion measurement found that the polar field through LaAlOs is 
much smaller than the theoretical predictions [42]. It is possible that various atomic recon- 
structions, such as cation intermixting, nonstoichiometry and defects could also significantly 
affect the electronic structures of these interface systems and be partially responsible for the 
discrepancies mentioned above. While a detailed discussion of atomic reconstructions is 
beyond our present work, it serves as an interesting direction for future research. 

Appendix A: Conduction electron and hole densities 

In this appendix, we explain how we calculate the conduction electron and hole densities 
for the symmetric superlattices and stoichiometric interface systems. The basic tool we use 
is the local density of states (LDOS) for the system 

D(t, E) = J2 \ip nk (r)\ 2 8(E - e nk ) (Al) 

rik 

which we integrate over the appropriate energy range. 
1. Double n-type and p-type superlattices 

The symmetric nonstoichiometric double n-type (p-type) superlattice is shown in Fig. 
[T^i (Fig. [lb). The corresponding electron (hole) denisty is presented in Fig. [2^ (Fig. [2b). 
For the symmetric superlattices, there is no polar field and analysis of the density of states 
shows that there is an easily identifiable energy gap at all locations in the film; the Fermi 
level is either above the conduction band edge for the the n-type superlattice or below the 
conduction band for the p-type superlattice. For the transferred electron density (Fig. 
we integrate the LDOS from the middle of the band gap to the Fermi level, 

D(r, E)dE (A2) 

where £ is the energy in the middle of the band gap. We are in fact counting all the electrons 
in the conduction bands. For the transferred hole density (Fig. [2b), similarly we integrate 
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the LDOS from the Fermi level to the middle of the band gap, 

e(r)= f D(v,E)dE (A3) 

i.e. counting all the holes in the valence bands. The integrated LDOS rj(r) and £(r) are 
then averaged over the xy plane and finally plotted along the z direction. 

2. Stoichiometric n-type and p-type interface systems 

The stoichiometric n-type (p-type) interface system is shown in Fig. [3b (Fig. |3b). The 
corresponding conduction electron (hole) density is shown in Fig. [2fc (Fig. |2H). For the 
stoichiometric n-type and p-type interfaces, the polar field ensures that the bands edges 
in the LaA103 are not flat and that the Fermi level will intersect the band edges of the 
SrTiC>3 as well as the surface of the LaA103 film (see Fig. |3]). However, the local density 
of states in the SrTiC>3 film still clearly shows a band gap. Therefore, we can also compute 
the transferred charges at the interface using the above formula. The only complication is 
that in addition to showing the transferred charges in the interface region, n(r) and £(r) 
will necessarily have contributions localized at the surface region of LaA103 which are not 
of direct interest when studying the interface alone; the choice of axis range in Fig. [2fc and 
EH effectively excludes this contribution. As a consistency check, we can also compute the 
transferred charge by computing the atomic projections of all Bloch states in the system 
and identifying all those states with the proper atomic character: Ti-ci or La-d character 
and partial occupancy for the n-type interface, and O-p character and partial occupancy for 
the p-type interface. One can then manually sum up these particular contributions to get 
the electron and hole distributions, and the results agree with the previous method in the 
relevant regions. 

Appendix B: Periodic boundary condition effects 

When we perform a slab calculation, the materials in the simulation are generally po- 
larized either due to intrinsic polar fields (e.g. LaAlOs) or external electric fields. In DFT 
calculations, the periodic boundary conditions (PBC) imposed on the simulation cell influ- 
ence the screening properties of the materials. We will show in the following two sections 
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that the presence of PBCs does not significantly affect the polarization due to intrinsic polar 
fields in the system, but it does induce large artificial effects when one studies the response of 
the system to an external field, requiring a careful set up of the simulation and interpretation 
of the results. 

1. Polarization due to intrinsic fields 

In this section, we illustrate that in the SrTiOs/LaA10 3 interface system, PBCs induce 
an artificial field through SrTiOs, but we show that this effect is neglegible. 

As shown in Fig. [T6a . the two boundary conditions of displacement D at each interface 
give: 

£>lao - £>sto = 4ttct (Bl) 

D STO -D y = (B2) 
PBCs require that the total potential drop over the whole simulation cell is equal to zero: 

E STO s + E LAO l + E Y v = (B3) 

where s, I and v are the thicknesses of SrTiOs, LaA10 3 and vacuum, respectively. Con- 
sidering that Ei = Di/ei {i is SrTi0 3 , LaA10 3 or vacuum), we can solve -D STO and .Dlao 
explicitly, 

D STO = -^a- j /£LAO (B4) 

£>lao = 47r(T— — ■ — — (B5) 

s/ e S TO + 1 1 clao + v/ e v 

The dielectric constants of SrTiOs, LaAlOs and vacuum are roughly: 

e ST o - 300 > e LAO ~ 25 > e v = 1 (B6) 

and the thicknesses of SrTiOs, LaAlOs and vacuum in the simulation cell are on the same 
order of magnitude (within a factor of two), therefore we have the following simplifications: 
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E LAO z> — (B7) 
£lao 

Ana I I 

^sto — = E LAO — > I-E'stoI |-Elao| (B8) 

eLAOCSTO v ve S To 

E y ~-—- = --E LAO -> ~ -^lao^ (B9) 
clao u u 

Eq. (IB8p and Eq. (1B9I) show that though SrTi0 3 is artificially polarized due to the PBCs, 
the electric field through SrTiOs is much smaller than the intrinsic field through LaAlOs and 
the voltage built across the LaAlOs film mostly drops in the vacuum. Based on the above 
approximate dielectric constants, the field in SrTi0 3 is only ~ 0.3% of that in LaAlOs. In 
addition, we can further reduce the field in SrTi03 by increasing the thickness of the vacuum 
v; the above relations give a quantitative measure of the (small) error for any finite value of 
v. 



2. Polarization due to an applied external electric field 

We illustrate in this section that when we use a slab calculation and a sawtooth potential 
to simulate an external electric field through nonpolar materials {e.g. SrTiOs), PBCs will 
artificially undermine the screening and give rise to a significantly smaller dielectric constant. 

Fig. fT6b shows a schematic of how an external electric field E ext along the z direction 
is screened in a nonpolar material. There are two induced electric fields. The one in the 
material is the stardard depolarization field Ed ep and the other field, which we call the 
outside field E out , is in the vacuum with the opposite direction. E ext , E dep and E out are 
the magnitudes of each field; their directions are explicitly shown in Fig. [16b. We denote 
the size of the unit cell along the z direction by L and the thickness of the slab by d. The 
sawtooth potential automatically satisfies the periodic boundary condition (the reversed 
part is not shown in Fig. [16b). Therefore, the induced electric fields required to satisfy the 
periodic boundary condition are given by: 

E dep d = E out {L - d) (BIO) 
On the other hand, Gauss's law gives: 
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FIG. 16: Schematics of how the periodic boundary condition affects the screening of the electric 
fields, a) Polarization due to intrinsic fields. The shaded parts are SrTi03 and LaA103, respec- 
tively. The interface is n-type. The unshaded part denotes vacuum, b) Polarization due to an 
applied electric field. The shaded part is a general nonpolar material and the empty part is vacuum. 



4ttS = E dep + E out (Bll) 
where the surface charge density X is related to the polarization in the material by 

£ = n-P = P (B12) 
and the polarization P is related to the total electric field E tot by 

P = xEtot (B13) 
where x is the permitivity. From Fig. [T6b . it is easy to see that 

E to t = E ext — E dep (B14) 

Combining Eq. (1B10HB14]) gives 



E ex t 



- 1 + 4, x( l-f) <B15) 



Now we can identify the dielectric constant e 
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t,rl l + 4irx(l-T) (BIG) 



E to t V 

Eq. ( IB 16 1) is different from the familiar formula e = 1 + 47r% in that we have an extra 
factor (1 — d/L). In typical slab calculations (1 — d/L) ~ 30% and for high-fc materials 
{e.g. SrTiOa) x ^> 1 so that e is dominated by the second term in Eq. ( 1B16I) and thus the 
reduction factor is a significant error. 



3. A remedy 

The origin of the deviation of the dielectric constant from the correct value of 1 + Attx 
is that the electric field induced by the bound charge ±S does not completely serve as the 
depolarization field, but is instead split into two parts, E dep and E out . The presence of E out 
is purely due to the imposed periodic boundary condition, and is unphysical. However, in 
practical slab calculations it is too computationally expensive to make d/L small. 

In order to get rid of E out , we need a simulation cell in which the depolarization field 
automatically satisfies periodic boundary conditions. Therefore, we use a mirror-symmetric 
simulation cell in which the external electric field is also mirror-symmetrically distributed. 
By symmetry, E&, p automatically satisfies periodic boundary conditions and E out is guar- 
anteed to be zero in the vacuum. Thus, the external electric field is correctly screened, but 
with a relatively modest increase in computational expense due to doubling the simulation 
cell along one direction. 



Appendix C: Field dependence of dielectric constant 



The Landau theory phenomenologically describes the field dependence of the dielectric 



constant 
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48| . It is assumed that the free energy of the system F(P, T) can be expanded 



in even powers of the polarization P: 



F(P, T) = F + AP 2 + BP 4 + CP 6 + ... (CI) 

where the coefficients A,B,C... may depend on the temperature T. Keeping only terms in 
F to the fourth order, we obtain: 
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and the permitivity follows: 



E = — = 2AP + 4PP 3 (C2) 
oP 



1 F)l? 

- = ^ = 2A + Y1BP 2 (C3) 
X dP 



Eq. ( 1C2j) and Eq. ( 1C3|) uniquely determine x — x{P) an d -P = P{E). The analytical solution 



to Eq. (1C2I) P = P(P) is complicated, but we can find a useful interpolation scheme 471. |48| 
by noting that as P is small, x ~ >* (2A) _1 and as P is large, x ^ -P 2 an d P oc P 3 , that is 

X oc P~ 2//3 , so 





(C4) 

For high-fc materials, we can also approximate e ~ Attx ^> 1, giving the final expression: 

(C5) 

\ V°o/ J 
where €q and are fitting parameters. 

We note that the truncation of Eq. flClj) is based on the assumption that P is small (i.e. 
E is small). However, the Landau theory itself does not give a characteristic polarization 
or electric field. Instead, we consider £ as a characteristic electric field whose value is 
determined by the fitting. Thus, data points much larger than So should not be used in the 
fitting because Eq. ( 1C1I) would break down. This is a self-consistent check. 



1. LaA10 3 

Before we calculate the dielectric constant of LaA103, we need to elucidate a subtle 
point. The we calculate here is only well defined in the thin slab of LaAlOs, not in the 
bulk. This is because we are in the region of very strong electric field (~ 0.7 V/olao)- With 
a few unit cells, the potential difference built across the slab will be larger than the band 
gap of LaAlOs, and Zener tunneling will occur. Bulk LaA103 can not accommodate such 
a large electric field and the Berry phase method |49] for calculating the dielectric constant 
breaks down in this regime. Instead, we resort to a mirror-symmetric slab calculation (see 
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• • • * • 

Sr Ti La Al O 



FIG. 17: Schematics of the simulation cell of LaA103 slab calculations. Both the simulation cell 
and the external electric field are mirror-symmetric. All fields through LaA103 are along the z 
direction (perpendicular to the interface). 

Fig. [T7]and Appendix |B]) . We turn on a mirror-symmetric external electric field E ext and 
calculate the macro-averaged electric field E to t through the material. The dielectric constant 
is then defined as: 



J^tot 

The raw data and the fitting curve using Eq. (1C5I) are shown in Fig. [181 The fitting results 
are: 



2. SrTi0 3 

Since SrTiC>3 is a nonpolar material with a large dielectric constant, the typical internal 
electric field through SrTiC>3 is much smaller than that through LaA103. The dielectric 
constant we are interested in can be defined in the bulk (based on the argument of metastable 
states) and calculated using the Berry phase method. In order to accurately determine 
the atom positions, we use 6 x 6 x 20 Appoint sampling and lower the force convergence 
threshold to 8 meV/ A. We directly calculate the total polarization P tot (both ionic and 
electronic) in the unit cell at a given total electric field E tot and the dielectric constant 
follows straightforwardly: 




(C6) 



e = 40.95, £ = 0.15V/ A 



(C7) 
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FIG. 18: Electric field dependence of dielectric constants of LaA10"3 and SrTi03. The solid squares 
are the results of the slab calculations and the solid triangles are the results of Berry phase calcu- 
lations. The solid lines are the fitting results using Eq. (|C5p . 



es(E } 



tot) 



. Ptot 1 



SIR 



(Q 



tot 



where fl is the volume of unit cell. The raw data and fitting curve are shown in Fig. ITS 
with the fitting results: 



4 = 309.6, £ 5 = 49.2V//im = 4.92 x 1(T 3 V/A 



(C9) 
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